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Nonionic contributions to the electric-field gradient at 181 Ta and m Cd impurity sites 
in R 2 3 (R= Sc, In, Lu, Yb, Tm, Er, Y, Ho, Dy, Gd, Eu, Sm) bixbyites 
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The time-differential perturbed-angular-correlation (TDPAC) technique was applied to the study 
of the internal electric-field gradient (EFG) in Eu- and Ho-sesquioxides in their cubic bixbyite 
phases. The results, as well as previous characterizations of the EFG at 181 Ta sites in oxides with 
the bixbyite structure, were compared to those obtained in experiments using 111 Cd as probe, and 
to point-charge model and ab initio results calculations for the EFG tensor at impurity sites in 
binary oxides. These studies provide quantitative information about electronic processes and the 
structural relaxations induced by the presence of impurity probes in the host lattices, and confirm 
the existence of nonionic contributions to the EFG in these systems. Our FP-LAPW calculations 
show that this nonionic contribution to the EFG is the dominating one, and that it is originated in 
the population of p states (5p in the case of Cd, 6p for Ta). 

PACS numbers: 61.72.Ww, 71.55.Ht, 81.05.Hd, 82.80.Ej 



I. INTRODUCTION 

During the last three decades, time-differential 7 — 
7 perturbed-angular-correlation (TDPAC) spectroscopy 
has been increasingly applied to condensed matter prob- 
lems through the precise characterization of the electric- 
field-gradient (EFG) tensor at diluted (ppm) radioac- 
tive probe atoms, adequately introduced in substitu- 
tional host lattice sites (see, e.g., Refs. 1-3 and refer- 
ences therein) . The characteristics of TDPAC allow car- 
rying out the experiments under a wide range of suit- 
able external conditions such as variable temperatures, 
pressures, atmospheres, etc. The single-atom counting of 
this technique, in combination with its highly localized 
sensitivity (due to the r~ 3 dependence of the electric- 
quadrupole interaction) , allows a detailed investigation of 
both structural and electronic properties of the systems, 
providing information about crystal structures, 4 crys- 
tal chemistry, 5 defects, 6,7 nanoscopic characterizations 
of surfaces, interfaces, 8 and highly dispersed species 9,10 , 
among other properties in solid state physics 11-18 , chem- 
istry and biology (see, e.g., Ref. 19, and references 
therein) . 

The very well suited ( m In — >) m Cd isotope is the 
most frequently used tracer in TDPAC experiments and 
has been largely applied to study semiconductor physics. 
In particular, after the initial work of Pasquevich et 
al. 20 on the internal oxidation of diluted indium impu- 
rities in Ag, a large amount of experimental work has 
focused on the EFG characterization at 111 Cd impurity 
sites in semiconductor and insulating binary oxides. 21 All 
the information that the EFG tensor can provide about 
the system under study could be obtained by confronta- 
tion of the experiment with an accurate prediction of 



the EFG, such as those obtained with ab initio calcu- 
lations. In the absence of such predictions at impurity 
sites, several attempts to correlate experimental results 
and semiempirical calculations have been made from the 
very beginning, 22,23 in order to describe the different con- 
tributions to the EFG at impurity sites. 24 In 1992, a well- 
defined empirical correlation between local (Cd valence 
states) and ionic contributions to the EFG was presented 
for nl Cd in binary oxides. 25 In that analysis, the local 
component of the EFG in binary oxides was extracted 
from all quadrupole-coupling constants measured at that 
time in TDPAC experiments with the m Cd probe. The 
resulting systematics revealed a linear dependence be- 
tween the local and ionic contributions to the EFG over 
a wide range of ionic EFG values. In 1994, Weht et 
al. 26 confirmed all the features of the correlation through 
an independent cluster calculation of the valence EFG 
contribution using the Extended Hiieckel Method. Re- 
cently, Errico et al. 27 ' 2a reported the first ab initio full- 
potential linearized-augmented-plane-wave (FP-LAPW) 
calculations of the EFG at an impurity site (Cd) in an ox- 
ide (Ti02), in excellent agreement with the experiments, 
supporting the existence of a dominant valence contribu- 
tion to the EFG in this system, which can be identified 
with the local contribution of Ref. 25. 

To study the influence of the electronic configuration 
of the impurity probe atom itself on the EFG it is es- 
sential to perform TDPAC measurements with different 
probes in isomorphus crystal structures. Examples of 
this kind of study are the work of Adams and Catchcn, 5 
and Shitu et al. 29 In TDPAC experiments, the second 
most commonly used radioactive probe is 181 Hf, which 
decays by /3~ to the 181 Ta isotope. Therefore, in 1990 
we started a comparative study of the EFG at 181 Ta and 
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m Cd sites in binary oxides with the aim of investigat- 
ing the EFG dependence on the electronic configuration 
of the probe atom and its coordination "chemistry" with 
its nearest neighbors. 30 Among the binary oxides, the 
sesquioxides that crystallize in the cubic bixbyite struc- 
ture constitute a group widely studied with the 111 Cd 
probe. 11 ' 12 ' 31-35 This group presents interesting proper- 
ties since it includes a large series of rare-earth as well 
as closed-shell metallic cations, a wide range of lattice 
constants (from 0.94 to 1.09 nm), and two nonequivalent 
cation sites with quite different symmetry, cation-oxygen 
bond lengths, and crystallographic abundance. These 
characteristics make this group ideally suited to study in 
detail the EFG dependence on the coordination geometry 
of the probe atoms and on their electronic configuration. 
With this purpose, in 1994 we started a systematic TD- 
PAC study of this group using 181 Ta as probe. 36-41 From 
this systematic study and using the same hypothesis as 
in Ref. 25, we found an empirical correlation between the 
local and ionic contributions to the EFG similar to those 
found for 111 Cd, but with different slope. 39 

In this work we report TDPAC experiments using 
181 Ta as probe in EU2O3 and H02O3 in order to com- 
plete the systematic study of the EFG at 181 Ta atoms 
located at defect-free cation sites in the bixbyite struc- 
ture. This systematics is compared to those previously 
established for the 111 Cd probe and to recent and new 
ab initio band-structure calculations of the EFG at im- 
purity sites in a small set of binary oxides. This compar- 
ison provides quantitative information about electronic 
processes and the structural relaxations induced by the 
presence of impurity probes in the host lattices, informa- 
tion that cannot be obtained (or is crudely estimated) by 
simple models such as the point-charge model, the use of 
antishielding factors, and arbitrary suppositions. The ab 
initio calculations also predict the existence of a dominat- 
ing nonionic (or local) contribution to the EFG originated 
from Cd-5p and Ta-6p states, in good agreement with the 
semiempirical model proposed in references 25 and 39. 
These new theoretical results, in combination with the 
experimental ones, enable us to discuss the validity of 
the widely used ionic model and to validate or discard 
the hypothesis used in the construction of the semiem- 
pirical model mentioned above and its predictions. 



II. SAMPLE PREPARATION AND TDPAC 
MEASUREMENTS 

Under suitable conditions, Fe, Mn, Sc, In, Tl, Y, and 
all the rare-earth elements form a sesquioxide. Polymor- 
phism is common among the rare-earth oxides and be- 
low about 2300 K three polymorphous crystallographic 
structures have been found: 42 the hexagonal A-, the mon- 
oclinic B-, and the cubic C-form (bixbyite). In the cubic 
structure the cations form a nearly cubic face-centered 
lattice (space group Ia3) in which six out of the eight 
tetrahedral sites are occupied by oxygen atoms. The el- 
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FIG. 1: Nearest-neighbor (NN) oxygen distribution around 
each cationic site (black atom) in the bixbyite structure. 



ementary cell of the oxide lattice consists of eight such 
cubes, containing 32 cations and 48 Oxygen ions. Two 
nonequivalent cationic sites, called C and D, both 06- 
coordinated, characterize the structure. Their relative 
abundance in the lattice is (/c//d) = 3. Site D is axi- 
ally symmetric and can be locally described (see Fig. 1) 
as a cation surrounded by six oxygen atoms localized at 
the corners of a distorted cube, leaving two corners of a 
diagonal free (T)3d point-group symmetry). At site C the 
cube is more distorted (C2-symmetry) and the six oxy- 
gen atoms leave two corners on a face diagonal free (see 
Fig. 1). 

In order to prepare the TDPAC samples, commer- 
cially obtained high-purity H02O3 and EU2O3 powders 
(Aldrich Products, 99.9% and 99.95% metallic purity, re- 
spectively) were treated in air for 24 h at 1023 K in or- 
der to achieve the crystalline C-phase and then pressed 
as circular pellets under 20 kN/m 2 pressure. Afterward, 
the H02O3 pellet was sintered in air for 2.5 h at 1273 
K whereas EU2O3 was sintered under similar conditions 
but at 1123 K for 1.5 h because in this oxide the C- 
form presents a phase transition to the monoclinic B- 
phase at around 1353 K. 42 In both cases, powder XRD 
analyses of the samples performed before and after the 
sintering showed that only the C-phase was present and 
also revealed an increase in the crystallinity of the sam- 
ples. The ion accelerator of the ISKP, Bonn, was then 
used to implant 181 Hf + ions into the samples with the 
following energies and doses: H02O3, 160 keV, 5xl0 13 
ions/cm 2 and EU2O3, 155 keV, 5xl0 12 ions/cm 2 . The 
as-implanted samples were subjected to TDPAC mea- 
surements in air at room temperature (RT = 300 K) 
and at atmospheric pressure. After these measurements 
the H02O3 and EU2O3 samples underwent an annealing 
treatment in air at 1323 K and 1123 K for 2 h, respec- 
tively, since in other rare-earth sesquioxides annealing 
above 1073 K for a couple of hours has demonstrated to 
be sufficient to restore the host crystallinity and to locate 
the impurities at defect-free cationic sites. In the case of 
H02O3 the TDPAC measurements were then carried out 
in air in the temperature range RT-1373 K in 100 K steps. 
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FIG. 2: R(t) spectra (left) and their corresponding Fourier 
transforms (right) of ( 181 Hf— >) 181 Ta in H02O3 measured at 
RT in the as-implanted stage and at the indicated tempera- 
ture after thermal annealing in air at 1323 K for 2 h. 



In the case of E112O3 the measurements were performed 
between RT and 1273 K (in 100 K steps) and at 77 K. 

The TDPAC technique is based on the determination 
of the influence of extranuclear fields on the correlation 
between the emission directions of two successive radi- 
ations emitted during a nuclear-decay cascade. A com- 
plete description of this technique can be found in the 
literature (see, e.g., Refs. 19 and 43). In order to perform 
the experiments presented in this work we made use of 
the well-known 133-482 keV 7 — 7 cascade in 181 Ta, pro- 
duced after the f3~ nuclear decay of the 181 Hf isotope. 
The TDPAC experiments were made using four BaF 2 
detectors in a coplanar 90° arrangement and a fast-fast 
logic coincidence system. The time and energy resolu- 
tions of the spectrometer for 181 Ta 7-rays are 0.65 ns 
(full width at half-maximum of the prompt peak) and 
9.5% for 7-rays of 662 keV ( 137 Cs), respectively. The 
experimental perturbation functions, R(t), were derived 
from eight concurrently measured coincidence spectra, 
four taken between detectors positioned with 180° sym- 
metry and four of eight possible with 90° symmetry. 37 To 
analyze the measured perturbation functions, a multiple- 
site model for nuclear-electric-quadrupole interactions for 
polycrystalline samples and spin 1=5/2 intermediate level 
of 181 Ta was used 
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FIG. 3: R(t) spectra (left) and their corresponding Fourier 
transforms (right) of ( 181 Hf^) 181 Ta in EU2O3 measured at 
RT in the as-implanted stage and at the indicated tempera- 
ture after thermal annealing in air at 1123 K for 2 h. 



R{t)=AH P G 22 {t) = 
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where fi is the relative fraction of nuclei that experiences 
a given perturbation and is the effective anisotropy 
of the 7 — 7 cascade. The u> n interaction frequencies are 
related to the quadrupole frequency luq — eQVzz/^Ofo 
by uj n ~ g n ^Q- The g n and S^n coefficients are known 
functions 44 of the asymmetry parameter 77 = (Vxx — 
Vyy)/Vzz, where Vu are the principal components of 
the EFG tensor that are arbitrarily labeled according to 
I Vxx I < \Vyy\ < \Vzz\- The exponential functions in 
Eq. 1 account for a Lorentzian frequency distribution of 
relative width S around ui n . 



III. 181 TA TDPAC RESULTS IN H0 2 3 AND 
EU2O3 

Figures 2 and 3 show the R(t) spectra and their corre- 
sponding Fourier transforms, taken on both as- implanted 
samples at RT in air, and the selected spectra taken after 
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annealing of the samples. Solid lines in the R(t) spectra 
are the best least-squares fits of Eq. 1 to the experimen- 
tal data. Solid lines in the Fourier spectra come from the 
Fourier transform of the R(t) fits. As can be seen from 
both figures, the as-implanted R(t) spectra are rather 
damped, as usually occurs in the as-implanted stage of 
181 Hf-implanted binary oxides, 36 ~ 38 - 40 showing that the 
radiation damage produces a not negligible host disorder 
around the probe atoms. After the annealing treatments, 
the radiation damage was removed in both samples. In 
effect, only two very well-defined (5 < 1.5%) hypcrfinc 
interactions (labeled Hfi C and Hfi D) were found, with 
the asymmetry parameter of Hfi D equal to zero and a 
high r\ value for Hfi C (see Table I) as predicted by the 
coordination symmetries of sites D and C, respectively. 
In addition, the frequency of the axially symmetric inter- 
action is twice as large as the asymmetric one, as found 
previously with the 181 Ta probe in Yb-, Y-, Dy-, Sm- 
, Lu-, Gd-, Er- and Tm-sesquioxidcs. 36-41 The asym- 
metry parameters and distributions of Hfi C and Hfi D 
for both oxides are constant with temperature, while the 
quadrupole frequencies show a reversible continuous de- 
crease with increasing temperature in the measured tem- 
perature range. The monotonic trend of the hypcrfinc 
parameters of Hfi C and Hfi I? is a typical behavior fol- 
lowed by hyperfine probes in substitutional cation sites 
of a single-phase crystalline structure. The total fraction 
of both interactions amounts to 100% of the probes in 
the whole temperature range of measurement for EU2O3. 
This statement also applies to H02O3, with the exception 
that two minor interactions amount to less than 20% of 
the probes in the temperature range RT-472 K. The rela- 
tive population of the two sites should be {fc/ Id) = 3 if 
the cationic sites were occupied according to their natu- 
ral abundance in the crystalline structure. In the present 
experiments, (fc/fo) = 3 was found for EU2O3. In the 
case of H02O3, the experimental ratio /c//d was smaller 
than this value. This departure was already observed 
for ( 181 Hf-^) 181 Ta in other bixbyites and was explained 
in terms of the relative ionic size of the probe and the 
cationic "space" of the host. 37,38 The present results for 
Ho- and Eu-sesquioxides are in excellent agreement with 
the occupancy trend previously observed (see Fig. 4). It 
is then clear that Hfi C and Hfi D can be undoubtedly as- 
signed to probes located at the crystalline nonequivalent 
cationic sites C and D. 

We want to mention that the /c//d ratio for H02O3 
and EU2O3 obtained with 181 Ta, as well as those for the 
rest of the bixbyite series, is independent of the sample 
temperature (in the range RT-1300 K). This may be ev- 
idence of no further ionic transport after the annealing 
process. 

The two minor frequencies observed for H02O3 exist in 
a short temperature range and disappear above 573 K in 
a reversible way. Due to the temperature range in which 
these minor interactions exist we have investigated the 
possibility of desorption-absorption of water in the sam- 
ples. However, a Differential-Thermal Analysis (DTA) 
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FIG. 4: Experimental ratio fc/fo as a function of the lattice 
parameter a for 181 Ta in bixbyites. 



carried out for the H02O3 samples did not confirm this 
hypothesis. 



IV. DISCUSSION 

A. Systematics of the EFG at 111 Cd and 181 Ta 
impurity sites in bixbyites 

In 1994, we started a TDPAC study of the EFG depen- 
dence on the probe atom in bixbyites, measuring the EFG 
at 181 Ta impurity sites in Yb-, Y-, and Dy-sesquioxides, 36 
ln 2 3 , 37 Sc 2 3 'and Sm 2 3 , 38 ' 39 Lu 2 3 , 39 Er 2 3 and 
Gd 2 3 , 4 ° and Tm 2 3 . 41 Including the new data pre- 
sented in this paper, i.e., Eu 2 3 and Ho 2 3 , we can 
now compare a set of twelve oxides, all having the same 
crystalline structure, but differing in their lattice con- 
stant by 10%. In Fig. 5 we present the values of ojq and 
■q for both sites as a function of the lattice parameter 
a. All the experimental data in this figure (and also in 
Fig 4) correspond to TDPAC experiments in which the 
181 Hf isotope was introduced into the samples by means 
of ion implantation, with the exception of Tni203. In 
this case the activity was introduced by means of a solid- 
state reaction assisted by ball milling between powders of 
Tm-sesquioxide and neutron-activated Hf0 2 . In spite of 
the limited accuracy of the derived hyperfine parameters 
in this experiment, due to a remaining 40% fraction of 
probes located in a highly distorted Hf02 phase, the hy- 
perfine interaction parameters corresponding to Tm203 
appear to be reliable considering their excellent agree- 
ment with the rest of the systematics. As can be seen 
in Fig. 5, there is a jump in the 181 Ta systematics. In 
effect, ujqc and uuqd for SC2O3 and In 2 03 (the sesquiox- 
ides with the smallest lattice parameters) depart from the 
linear dependence on a shown by the rest of the sesquiox- 
ides. Concerning the asymmetry parameter, while r\u is 
nearly independent of a for all the measured sesquiox- 
ides and close to 0, r\c does not show an increase with 
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TABLE I: Results of least-squares fits of Eq. 1 to the R(t) spectra displayed in Figs. 2 and 3, taken at RT after annealing. /, 
and Si are expressed in % and ljq; in Mrad/s. When no errors are quoted it means that the parameter was kept fixed in order 
to evaluate errors. In the case of H02O3 the fractions fc + /d are normalized to 100% (see text). 
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FIG. 5: Experimental quadrupole frequency loq and asym- 
metry parameter 77 for substitutional 1 Ta impurities at sites 
C (black) and D (hollow) in bixbyites plotted as a function 
of the lattice parameter a at RT. To obtain V zz from loq 
we used Qr a = (+)2.36(5) 6. 45 The errors are smaller than 
the symbols. Dotted lines are the PCM predictions. The loq 
predictions are normalized to the In2C>3 values. The lattice pa- 
rameters and the atomic coordinates used in the PCM calcu- 
lations for each oxide have been taken from neutron and x-ray 
diffraction (XRD) determinations quoted in the literature. 46 
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FIG. 6: Experimental quadrupole frequency loq and asymme- 
try parameter r\ for substitutional 111 Cd impurities at sites C 
(black) and D (hollow) in bixbyites 11 ' 12 ' 31 " 35 plotted as a 
function of the lattice parameter a at RT. The results for 
LU2O3 and Tm203 correspond to 1273 and 583 K, respec- 
tively (see Refs. 11 and 35). In order to obtain V zz from loq 
we used Qcd = (+)0. 83(13) b . 47 The errors are smaller than 
the symbols. Dotted lines are the PCM predictions. The loq 
predictions are normalized to the In2 03 values. 



a for In- and Sc-sesquioxides, as occurs for the rest of 
the systematics (see Fig. 5). The anomalous behavior 
of the loq and r\c parameters may indicate the presence 
of structural distortions in the host lattice induced by 
the impurity. We will discuss this statement in the next 
section. 

If we compare the trend of loq and r\ with point-charge 
model (PCM) predictions (see Fig. 5), we find that the 
observed increase of r\c and the dependence of loqc and 
loqd on a are in poor agreement. At this point it is nec- 
essary to analyze the experimental hyperfine interaction 
parameters at 111 Cd impurity sites in bixbyites reported 
in the literature. 11 ' 12,31 ~ 35 Figure 6 shows a monotonic 



dependence of the hyperfine parameters loq and 77 on the 
lattice parameter a for both sites. These trends are in 
poor agreement with those predicted by the PCM, as 
can be seen in the same figure. In the past, this discrep- 
ancy was explained by two different approaches. In the 
first one, 33,48 the PCM was considered accurate enough 
to exactly predict the hyperfine parameters that charac- 
terize the EFG at Cd sites in bixbyites. In this model, 
the EFG is given by 

V^ ic = (l-7oo)V&* ( 2 ) 

where 700 is the Sternheimer antishielding factor and, in 
the principal axes system: 
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In the above expression, Zk denotes the (arbitrary) ionic 
charge, Xki and Xkj the coordinates, and the distance 
of the k th ion to the site where the EFG is calculated, 
located at the origin of the coordinated system for sim- 
plicity. The discrepancy between PCM predictions and 
the experimental results was attributed to a wrong de- 
termination of the atomic positions in these compounds 
(mainly the oxygen atoms). Within this framework, the 
authors of Refs. 33 and 48 used the PCM to refine the 
atomic positions in order to reproduce the experimental 
EFG results. It is worthy of mention that that approach 
presented several assumptions, e.g., pure ionic bonding 
between the probe atom and its neighbors was assumed 
and local lattice distortions caused by the 111 Cd impu- 
rity were not considered. At this point, it is important 
to mention that that model was proposed before a 181 Ta 
systematics in bixbyites was established. Now, the very 
different trend between the experimental 181 Ta data pre- 
sented here and the PCM predictions with the refined 
coordinates of Ref. 48 gives clear evidence of the failure 
of this model in the 181 Ta case and makes its applica- 
bility suspicious in the m Cd one. In the next section, 
we present theoretical evidence that supports this state- 
ment. 

The second approach was first proposed to describe 
the whole set of EFG results at 111 Cd in binary oxides in 
a unified way, 25 and then it was extended to the 181 Ta 
impurity. Since this model was constructed to exactly 
reproduce the magnitude of all the existent experimental 
EFG results, the hypothesis involved in it and the predic- 
tions of this model for the sign of V zz (opposite to those 
of V'z Z %c , see Ref. 25) have to be discussed now, enlight- 
ened by the first- principles predictions, which were not 
available at the time the model was proposed. In this 
approach, Rcnteria et al. 24 ' 25 proposed the existence of a 
dominating local contribution to the EFG. In this model 
the EFG is given by 



Vzz = V l z T ' + (l-7co)T4z • 



(4) 



The local contribution to the EFG takes into account 
covalence between the probe and its neighbors. From 
the study of the experimental EFG results at 111 Cd and 
181 Ta impurity sites in binary oxides we demonstrated 
that, for both probes, V^£ al shows an almost linear de- 
pendence on V^z %c ( see R- e f- 39). We also demonstrated 
that Vz°z al can be factorized in two parts, one that de- 
pends on the geometry of the cationic site and another 
that depends on the electronic configuration of the probe, 
thus 39 - 40 



Vzz (site i, probe i) 



probe i 



(1 - 7oo) 



probe i 



X geometry y 1 ^ {site i). (5) 



We must say here that this semiempirical model is 
still based on point-charge summations and takes into 
account the contribution of the electronic structure of 
the probe to the EFG by means of antishiclding factors 
and empirical parameters (such as ^ rohe and ^9 eometr vy 
Moreover, this model does not take into account possi- 
ble effects introduced by the impurity character of the 
probe, such us lattice distortions. Nevertheless, this sim- 
ple semiempirical model has been successfully used as a 
tool for the assignment of hypcrfine interactions in com- 
plex systems. 16 In the next section the hypothesis and 
predictions of this model arc discussed in terms of very 
recent band-structure FP-LAPW calculations of the EFG 
at Cd and Ta impurity sites in In 2 03 and other binary 
oxides. 



B. FP-LAPW results for the EFG tensor at 
impurity sites in binary oxides 

Since this paper deals with the EFG at impurity sites, 
the interpretation of the experimental results involves 
the understanding of chemical differences between the 
probe atom and the indigenous ion replaced by the im- 
purity. The experimental results show that the differ- 
ences between probes and indigenous atoms are manifest 
in subtle ways that are not well described by conven- 
tional models, as already pointed out in Ref. 5. For 
an accurate calculation of the EFG, the electronic con- 
figuration of the host, perturbed by the presence of the 
impurity, has to be determined. This can be done in the 
framework of the Density- Functional Theory (DFT). In 
this kind of calculation, electronic and structural effects 
introduced in the host by the presence of the impurity 
probe (impurity levels, structural distortions, etc.) can 
be described without the use of arbitrary suppositions. 
Unfortunately, these calculations prove to be not trivial 
and time-consuming. For this reason, very few calcu- 
lations have been performed in systems with impurities 
and the method is far from being routinely applied in 
this field. 27 ' 28 Moreover, the applicability of the DFT to 
bixbyites that contain rare earth cations is not clear yet 
since the 4/-valence electrons are not well described by 
this theory. 49 

In the last years we performed FP-LAPW 50 cal- 
culations in the systems Ti0 2 (Cd), 27 < 28 Sn0 2 (Cd), 51 
Ti0 2 (Ta), and Sn0 2 (Ta). 52 We present here new calcu- 
lations of the EFG tensor at Cd and Ta impurities that 
replace cations in the model case In 2 03. For the cal- 
culations we employed the Wien97.10 implementation 53 
of the FP-LAPW method within the Generalized Gradi- 
ent Approximation (GGA) for the exchange-correlation 
potential. 54 For the parameter RKmax, which controls 
the size of the basis set, we took the value of 7 and in- 
troduced local orbitals for In-4d and Ap, 0-2s, Cd-4rf 
and 4p, and Ta-6s, bp, and 4/ orbitals. With these pa- 
rameters, the convergence errors in the Vu components 
are smaller than 0.1 x 10 21 V/m 2 . In the case of the 
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asymmetry parameter 77, the convergence error is smaller 
than 0.05. In order to study the structural distortions 
introduced by the impurity in the host we considered the 
displacements of the impurity nearest oxygen neighbors 
until forces acting on these atoms vanished, assuming 
that structural distortions preserve the point-group sym- 
metry of the cell in its initial configuration. Due to the 
large number of atoms in the unit cell of In 2 03, we did 
not use a super-cell approach for the simulation of the 
isolated impurity. Details of the method of calculation 
and the way to deal with the impurity were extensively 
explained in Refs. 27 and 28. It is clear that, even includ- 
ing the new FP-LAPW results reported here, the number 
of EFG calculations at impurity sites is very small to per- 
form a systematic study, but some important conclusions 
can be drawn. 

The first important theoretical result is that in all the 
studied systems the impurity probes introduce local dis- 
tortions in the host. It is worthy of mention that in all the 
cases the amount of distortion per atom decreases rapidly 
from the impurity nearest oxygen neighbors (Onn) to 
other shells. In the case of Cd located at both cationic 
sites of In2C>3 we found that the impurity produced relax- 
ations of about 5% (see Table II) of the cadmium Onn- 
Additionally, in the case of Cd located at site C we found 
that the impurity moved out from the symmetry site, but 
this displacement (3xl0~ 4 nm) is very small compared 
to the displacement of the Cd Onn- Larger relaxations 
were also found in Ti02 and Sn02 for the Cd impurity. 
These relaxations can be understood from the fact that 
the bond length of the sixfold coordinated Cd ion in CdO 
is about 0.235 nm. It seems that the local structure of the 
Cd impurity tries to reconstruct the environment of Cd in 
its own oxide. 51,52 Since in the rest of the bixbyites most 
of the bond lengths are smaller than 0.235 nm, one could 
also expect relaxations around the Cd impurity site in 
all of them. The presence of these structural relaxations 
in bixbyites with Cd shows that the refinement proposed 
in Refs. 33 and 48 is not valid since it is based on the 
fact that the Cd impurity does not introduce structural 
distortions in the host. In this sense this model mistakes 
a local effect (the O nn relaxations originated by the im- 
purity character of the probe) for a pretended wrong de- 
termination of the host atomic positions by XRD, which 
involves all the lattice. 

In the case of Ta at site D of 1^03, we found that 
Ta induces a 7% contraction in the Ta,-0 nn distance. 
Calculations at site C predict 5-10% contractions of the 
Ta-0 nn distances (see Table II) and a negligible dis- 
placement of Ta from the symmetry site. We have also 
found contractions in the case of Sn02(Ta), whereas in 
Ti02(Ta) the local structural distortions are negligible. 52 
These contractions can be understood from the fact that 
the bond lengths of the sixfold coordinated Ta ion in 
Ta02 are about 0.202 nm. As in the case of Cd, it seems 
that the local structure tries to reconstruct the environ- 
ment of Ta in its oxide. Due to the larger bond lengths 
of the rest of the bixbyites we should also expect con- 



tractions of the Ta-OATAr distances. From these consider- 
ations it is clear that distortions with different symmetry 
and/or magnitude for In- and Sc-sesquioxidcs with re- 
spect to those of the rest of the systematics could be 
at the origin of the anomalous behavior of the hyperfine 
parameters for Ta in bixbyites shown in Fig. 5. 

We can discuss now the FP-LAPW predictions for the 
EFG tensor and compare them with the experiments, 
the PCM, and the semiempirical model of Ref. 39. From 
our calculations we found that relaxations (or contrac- 
tions) beyond the Onn shell do not produce qualitative 
changes in the EFG tensor. For this reason all the FP- 
LAPW results for the EFG tensor shown here correspond 
to the relaxed position of the Ojvjv. There is an excellent 
agreement between the experiments and FP-LAPW pre- 
dictions for the EFG at Cd impurity sites, as shown in 
Table III. On the other hand, the PCM fails completely 
in the description of the EFG for Ti02- The agreement 
is better for Sn02, but the PCM fails in the prediction 
of the asymmetry parameter. This failure should be ex- 
pected due to their covalency originated in the short Cd- 
O nn bond lengths. In the case of In 2 03, where the bond 
lengths are larger, the PCM agrees with the FP-LAPW 
EFG prediction for Cd at site D, but disagrees in the 
sign of the major component of the EFG tensor for Cd 
at site C. 

In the case of the Ta impurity, the agreement between 
the experiments and FP-LAPW predictions is very good, 
with exception of Ta at site C in 1^03. Here it is impor- 
tant to say that due to the hypothesis used in the relax- 
ation process (structural relaxations preserve the point- 
group symmetry of the cell in its initial configuration), 
FP-LAPW cannot reproduce the small departure from 
the axial symmetry experimentally observed for Ta at 
site D of tri203 (see Table III). On the other hand, the 
PCM fails completely in the description of the EFG at 
Ta impurities located at cationic sites of Ti02 and Sn02- 
In the case of the EFG at Ta impurities located at sites C 
and D of In 2 03, the PCM agreement with the available 
experimental data (Vzz and rf) is better than that of the 
FP-LAPW predictions. Hence, from eight probe-oxygen 
configurations, five are not well described by the PCM. 

The ab initio calculations predicted local distortions 
in almost all the cases studied. If we use the relaxed 
positions of the Oat n atoms predicted by FP-LAPW in 
the PCM calculations, this model still fails in the pre- 
dictions of the EFG. Moreover, in most of the cases the 
disagreement becomes worse. From this discussion it is 
clear now that the failure of the PCM prediction of the 
EFG at impurity sites is not only due to the use of wrong 
atomic positions but to a wrong description of the elec- 
tronic structure of the impurity-host system, in particu- 
lar in the vicinity of the impurity. 

Based on the FP-LAPW calculations we can finally dis- 
cuss the validity of the hypothesis on which the semiem- 
pirical model is based, i.e., the existence of a dominating 
local contribution to the EFG. In this model, the EFG 
can be written as 
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TABLE II: Experimental and FP-LAPW-refined In-O distances (in nm) in pure In 2 03 compared with the relaxed Cd-0 and 
contracted Ta-0 distances (last two rows) calculated when a Cd or a Ta atom replaces an In atom in ln 2 3 . 





Site D 
d(In - O) 


d(In - Ol) 


Site C 
d(In - 02) 


d(In - 03) 


Pure ln 2 3 (Experimental) 513 


0.219(2) 


0.212(2) 


0.219(1) 


0.221(2) 


Pure ln 2 3 (FP-LAPW, this work) 


0.218 


0.213 


0.219 


0.223 


In 2 3 :Cd (FP-LAPW, this work) 


0.227 


0.220 


0.229 


0.233 


In 2 3 :Ta (FP-LAPW, this work) 


0.202 


0.197 


0.197 


0.209 



TABLE III: Experimental EFG results available for Ti0 2 , Sn0 2 , and ln 2 3 doped with Cd and Ta compared with the PCM 
and FP-LAPW predictions. The signs of Vzz at the Cd and Ta sites are not known experimentally. 



System 


Method 


Vzz 


Vzz direction 


n 


Ti0 2 : Cd 


PCM 


—2.21 


[1101 


0.43 




FP-LAPW 28 


+4.55 


[1101 

L J 


0.26 




Experimental 27 


5.34(1) 


[110] or [110] 


0.18(1) 


Sn0 2 : Cd 


PCM 


+5.21 


[1101 


0.40 




FP-LAPW 51 


+5.35 


[1101 
L J 


0.10 




Experimental 56 


5.83(4) 




0.18(2) 


In 2 3 :Cd (Site D) 


PCM 


+7.5 


[111] 


0.00 




FP-LAPW (this work) 


+7.6 


[111] 


0.00 




Experimental 13 


7.7(1) 


[lll] a 


0.00(5) 


In 2 3 :Cd (Site C) 


PCM 


-4.9 


[00.751] 


0.83 




FP-LAPW (this work) 


+5.6 


[00.751] 


0.68 




Experimental 13 


5.9(1) 


[0Tl]° 


0.69(1) 


Ti0 2 :Ta 


PCM 


-4.55 


[110] 


0.43 




FP-LAPW 52 


-13.0 


[001] 


0.69 




Experimental 57 


13.3(1) 


[001] 


0.56(1) 


Sn0 2 :Ta 


PCM 


+10.76 


[110] 


0.40 




FP-LAPW 51 


-17.7 


[001] 


0.96 




Experimental 30 


17.04(1) 




0.70(2) 


In 2 3 :Ta (Site D) 


PCM 


+15.46 


[111] 


0.00 




FP-LAPW (this work) 


+19.7 


[111] 


0.00 




Experimental 37 


16.31(2) 




0.149(7) 


In 2 3 :Ta (Site C) 


PCM 


-10.10 


[00.751] 


0.83 




FP-LAPW (this work) 


-17.8 


[010.3] 


0.26 




Experimental 37 


12.35(9) 




0.60(9) 



a Experimental values obtained with 111 Cd in Er 2 3 and Ho 2 3 
single crystals. 33 



r latt{PCM) 



^Aatt(PCM) tc -, 

Vzz > (6) 



Vt /-Local t rtaz 

zz = V zz - 7oo V zz 

meanwhile in the FP-LAPW framework we have 



j/LAPW _ irval 
V ZZ — V ZZ 



V, 



latt(LAPW) 



ZZ 



-irval,p , Trval,d . -.Aatt(LAPW) 

— ZZ ' v 77. • " 7 



ZZ 



' ZZ 



(7) 



where V ZZ ' P (V ZZ ' ) is the valence contribution of Cd- or 
Ta-p(d) orbitals inside the muffin-tin sphere of the impu- 
rity and y l ^-^ LAPW ) ; s t ne total contribution originated 
in the charge density outside that sphere. 

In both cases (Eqs. 6 and 7) the lattice contributions 
are negligible. Taking into account the d orbitals involved 
in the 7oo factor (4c? for Cd, 5d for Ta), we can correlate 



the term -700 y l ^^ PCM ) f jt^ 6 with the d valence con- 
tribution to the EFG in FP-LAPW calculations (Eq. 7). 
In this way, the p valence contribution of Eq. 7 has to 
be correlated with V z °£ al in Eq. 6. With this "bridge" 
linking both models, we can return to the results ob- 
tained in the FP-LAPW calculations. In all the studied 
systems (including 1^03) and for both probes we found 
that the p contribution (5p in the case of Cd, 6p for Ta) 
dominates over the d contribution. Based on these re- 
sults, we can conclude that the ab initio calculations of 
the EFG tensor at impurity sites support the hypothesis 
of the existence of a local and dominating contribution 
to the EFG on which the semicmpirical model is based. 
This model gives a better qualitative description of the 
EFG than the PCM because it includes a contribution 
to the EFG that is not taken into account in the PCM. 
But the semiempirical model is based on the use of the 
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7oo factor and point-charge summations that do not take 
into account the structural distortions introduced by the 
impurity. These shortcomings could produce wrong pre- 
dictions in our model, such as the wrong prediction of 
the EFG sign in various of the compounds shown in Ta- 
ble III. An extensive first-principles study of the EFG 
at impurity sites in binary oxides is essential in order to 
remove these shortcomings. 

V. CONCLUSIONS 

In order to compare complete sets of EFGs at 181 Ta 
and 111 Cd sites in bixbyite oxides we determined the 
electric-quadrupole hypcrfine interactions at 181 Ta im- 
purities in EU2O3 and H02O3. Including previous re- 
sults, we could compare a set of twelve oxides having the 
same crystalline structure. Our study shows that in all 
these oxides, after the implantation and subsequent high- 
temperature annealing, the probes 181 Hf occupy two sites 
that have been identified as the two nonequivalent cation 
sites of the bixbyite structure. We confirm a jump in the 
values of luqc and luqd that takes place for a < 1.0391 
nm, that is for lattice parameters smaller than that of 
Lu-sesquioxide. Concerning the asymmetry parameter, 
while t]d 0.0 is nearly independent of a for all the mea- 
sured sesquioxides, r\c for indium and scandium sesquiox- 
ides departs from the tendency that rules its behavior 
for the rest of the systematics. These anomalous exper- 
imental behaviors, in combination with our FP-LAPW 
predictions for Ta in L12O3 (and in the other binary ox- 
ides), suggest that Ta induces contractions of its Onn 
in all the bixbyites, but with different symmetry and/or 
magnitude in the case of In203 and SC2O3. 

The experimental EFG data obtained at 111 Cd and 



181 Ta impurity sites in binary oxides available at present 
clearly show that a single ionic model cannot describe 
the EFG, suggesting the existence of a dominating non- 
ionic contribution as stated in the semiempirical model. 
This assertion is supported by our band-structure FP- 
LAPW calculations at Cd and Ta impurity sites in binary 
oxides, which show that the dominating contribution to 
the EFG is originated in Cd-5p and Ta-6p states, respec- 
tively. These results can now explain the dependence of 
the local EFG, proposed in the semiempirical model, on 
the electronic configuration of the probes. Future FP- 
LAPW calculations, currently in progress, at impurity 
sites in oxides and, in particular, in bixbyites, will cer- 
tainly go deeper into the understanding of structural and 
electronic properties at impurity sites in semiconducting 
oxides. 
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